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ABSTRACT 



The purpose of this thesis was to evaluate the accuracy of the existing Mini- 
Ranger network in Monterey Bay and to suggest means by which the network 
can be improved. This network consists of six stations located around the Bay 
and installed by the Monterey Bay Aquarium Institute (MBARl) and the Nat'al 
Postgraduate School (NPS). 

In order to undertake this task, data were made available by MBARl which 
were collected on a cruise made by their vessel "Pt. Lobos". Additional data were 
gathered on a second cruise by the \'essel "Pt.Sur" of the NPS. 

The data analysis indicated network problems. An effort was made to iden- 
tify these problems and to compute \arious correctors. In addition, an equation 
has been derited which enab’es use of Mini-Ranger data collected when signal 
strengths are low. 

The estimation of the accuracies obtained from the network through the t'ar- 
ious tests applied, and the conclusions drawn, can be used as a guide to future 
users of the svstem. 
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I. INTRODUCTION 



In 1988, the Monterey Bay Aquarium Research Institute (MBARI) and the 
Naval Postgraduate School (NPS) began work to create a permanent network of 
Mini-Ranger stations around Monterey Bay. The net was to provide precise po- 
sitioning for research vessels and other scientific operations. Initially MBARI 
was seeking coverage of as much of the bay as possible, especially the Monterey 
canyon feature, with absolute horizontal position acuracies on the order of a few 
meters. It was understood this would be difficult to achieve, given the equip- 
ment's specified range (40 nm). This study was undertaken to determine what 
levels of accuracy and spacial coverage were being achieved, as the net was con- 
figured in early 1989, and also to suggest possible improt ements. 

Positioning data from two cruises were examined, one from a "Point Lobos" 
cruise on 31 MAR 89. and a second from "Point Sur" cruise on 22 SEP 89 (the 
"Point Lobos" and "Point Sur" are vessels used by MBARI and NPS, respec- 
tively). Significant errors were found in a high percentage of the position fixes 
from the "point Lobos" data set. It was apparent that at least some of the prob- 
lem arose from the use of weak signals (low signal strength readings) and incom- 
pletely calibrated equipment. Consequently, additional measurements were made 
onshore o\er known length baselines in order to quantify the systematic and 
random error components of the range measurement under various signal 
strength conditions.! These findings stimulate the direction of the rest of the 
study: 

• What were the main sources of errors? 

• Is there any postprocessing that can correct the data and derive better posi- 
tions? 

• What is a realistic estimate of the accuracy achiex’ed? 

• Are there any chances for improvement? 

1 It should be noted that the "Point Lobos" .Master station failed before the baseline testing 
was run. 1 hus it was not possible to confirm the error estimates derived from analysis of the cruise 
data with independent baseline measurements. 
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II. DESCRIPTION 



A. MINI-RANGER SYSTEM 

The Motorola Mini-Ranger system being used by MBARI and NFS is the 
Falcon 484. It is widely used for precise positioning in a t'ariety of nearshore 
marine applications like hydrographic surveying, dredging, oceanographic data 
collection, etc. Numerous excellent descriptions of the system are available in the 
literature. 

The system operates on the principle of a pulse coded transponder. The 
Master station (also refered as to as the Recei\’er-Transmitter) is located on the 
ship and ranges to a maximum of four Reference stations ashore. The maximum 
update rate is once per second. The Master and Reference stations transmit on 
different frequencies in the microwa\e band. 5570 and 5480 Mhz, respecti\’ely. 
Both range and signal strength data are output for each Reference station inter- 
rogated. 

Signal strengths are computed at the Master as a relative indication of the 
power level recei\ed from the Reference station. Values can ^•ary from about 5 
to 99. The manufacturer warns that any ranges with signal strength below 13 are 
"weak" signals, subject to random errors exceeding the 2 meters (la) specifica- 
tion. Signal strength ob\ iously falls off with increasing range, but can also drop 
due to any degradation of propagation conditions between the ship and shore 
station. Hence, weak signals can and do occur at any range. Monitoring of the 
received signal strength is critical to assessing the accuracy of all range data. 

The range is computed at the Master station based on the measured round- 
trip travel time for the pulse and a preset refractive index for microwave propa- 
gation in the atmosphere. Unless changed by the user, the receiver assumes a 
value of N = (n-l )xl 0^ = 320 as an indicator for the refracti\ e index (n), which is 
the value used throughout this study. Output ranges also can be corrected for any 
constant (systematic) error previously entered by the user. Because this constant 
error includes time delays at both ends of the line, it is essential to ha\ c a specific 



calibration of each Master/Reference station pair. The manufacturer notes that 
without calibration, systematic range errors of up to ten meters are possible. 

It is possible to have the receiver correct ranges from "slope to horizontal" 
distances before output, based upon user provided elevations of the Master and 
Reference stations. During this study, both the MBARI and NFS equipment 
were configured to output "slope" ranges only, choosing to handle this correction 
as part of the position computation algorithm. 

The Reference stations are supplied with electric current of 22 to 32 Volts 
DC from a pair of common car batteries, or permanently from an AC supply 
through a converter. A possible insufficient current supply can cause errors in the 
displayed distance in the receiver [Mini-Ranger Operation Manual, p.II-1 1]. 

The standard system has a maximum range of 37 km (20 n.m). with Refer- 
ence station antennas of 13 db. The maximum range can be extented to 75 km 
(40 n.m) by raising the antenna gain of the Reference stations to 19 db, if the 
station's ele\ ation is sufficient (due to the dip of the horizon). The master station, 
which is on the \essel. has an omni directional antenna of 6 db. The discussion 
and results in this thesis are all given using e.xtended (19 db) gain antennas on the 
shore stations. 

Another significant factor for a Reference station is the directi\ ity of its an- 
tenna. The gi\ en maximum range of the system corresponds to the main lobe di- 
rection and drops as one mo\ es to the sides. Figure 1 shows the pattern of a 19 
db gain antenna of a shore station. We can see that the range drops under 30 
km at 40*’ off the center direction. This is significant for a station that uses its 
maximum range ability and thus makes the choise of antenna directi\ ity critical. 

B. MONTEREY BAY NETWORK DESCRIPTION 

The Motorola Mini-Ranger network in Monterey Bay consists of 6 stations 
with high gain antennas (19 db). installed permanently on known positions 
around the Bay. The installation was done by both MBARI and NFS serving 
primarily the Biological and Oceanographic research of their \essels "Foint 
Lobos" and "Foint Sur" respecti\ely. The positions at which the stations are in- 
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Figure 1. 





19-db ANTENNA PATTERNS. 
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Table I. PRELIMINARY POSITIONS OF THE MINI-RANGER 
STATIONS 



Name 


Latitude(N) 
Longitudc(W) 
deg min sec 


UTM coordi- 
nates 


Scale 

factor 


Ele- 

vation 

(m) 


Trevor 


36 58 28.785 
122 02 31.536 


585260.161 

4092490.284 


.999690 


51.7 


Wats 


36 54 26.895 
121 50 39.819 


602947.672 

4085231.399 


.999731 


23.1 


Packard 


36 49 44.407 
121 46 04.912 


609863.128 

4076611.345 


.999749 


33.9 


Hays 


36 38 33.826 
121 47 45.895 


607621.289 

4055915.264 


.999743 


137.2 


Doppler 


36 36 05.601 
121 52 37.479 


600434.009 

4051260.014 


.999724 


10.0 


Hank 


36 36 25.381 
121 55 08.722 


596669.490 

4051826.411 


.999715 


134.5 



Stalled are not ideal due both to the permanent power supply problem and to the 
fact that some of the ideal station positions are on pri\ ate properties. 

On board the "Point Lobos" the Mini- Ranger omni-directional antenna was 
installed on the mast 9.1 meters above the water level. Thus the maximum 
achie\ able range due to the line of sight effect was increased. 

The shores of the Bay and the nearshore areas are generally low sand dunes 
so some stations don't ha\'e the proper altitude to reach their maximum range. 

Table 1 gi\cs the preliminary positions [Schncbelc. 1989] of the stations, the 
UTM grid coordinates and their scale factor. The horizontal datum is the 
American Datum of 1983 (NAD 83). The UTM projection used is for Zone 10 
(central meridian 123'’ W.) and the scale factor 0.9996. Any errors in the posi- 
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lions of these stations will propagate into the computed vessel position by an 
amount that depends on the geometry. If the errors are smaller than several 
decimeters then there is no significant effect. Revised positions are shown in Ta- 
ble 3, section IV-C. 

The stations are not ideally positioned and so may be moved to different po- 
sitions in the future to allow better coverage of the areas of interest. 

Station Doppler is not yet permanently installed for safety reasons. 
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III. THEORETICAL APPROACH TO THE PROBLEM 



A. THE RANGING PROBLEM 

The technique used for computing the ship's position using ranges from 
known reference stations will be derived and explained in the beginning so the 
reader will be able to follow and understand the data manipulation and results 
of this study. It's well known that one way to handle the computations involving 
ranging data is to work on some appropriate projection (e.g the UTM projection). 
There are se\eral techniques that can be used for the position derivation: 

• Ateraging the coordinates resulting from all possible combinations of meas- 
ured ranges. 

• Position calculation using the two most accurate ranges (using of the rest as 

a check). 

• Graphical methods. 

In this thesis we use the t’ariation of coordinates method which is described 
in Cross [1981]. For the position derivation, at least 2 ranges are needed (two 
unknown parameters in the plane solution X.Y). With more ranges one can esti- 
mate the accuracy of the resulting position and also make statistical decisions 
about the quality of the data. 

Figure 2 illustrates the concept of ranging from the i-th station. Consider a 
station S. in a coordinate system XY. The obsert'ed range, reduced to the plane, 
is OR, (not shown). Consider an assumed position (AP). The distance between 
the reference station and the assumed position is CR,. The calculation of the 
values dX.dY. which represent the corrections that hate to be applied on the AP, 
will give the desired true position (TP) after some iterations. 

As is obvious from the figure: 

= ^AP ~ ''^Si 
( i)AP = ^ AP~ ^ Si 

The key equation is the one for the calculated range: 
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Figure 2. VARIATION OF COORDINAI ES METHOD 



- {x^Ap + (5 ()‘( 7 > 



In order to find the transformation from range coordinates to position coor- 
dinates it is necessary to find the derivati\e: 



2CR,dCRi = 2(A’,),,pc/A 



dCRi = MdX + NdY 



where 
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The estimate of the true range then is equal to the observed range plus an 
error vv^hich can be positive or negative, called the residual V,. It is also equal to 
the calculated range plus the increment ciCR,. 

The true range then is OR^+l] where F, is the residual. But it's also 
CR, + dCRi. Solv ing for dCR„ we take: 



For the n ranges problem, in matri.x form, the equation may be slated as fol- 
lows: 



dCRi = [OR - CR)i + 1) 



and finally; 



M,dX + N,dV = [OR - CR)i + f',- 



AX= B + V 



Where matrices A ,B and X respectively are: 







n 



N. 



n 



{OR - CR)^ 
{OR-CR)^ 

{OR - CR), 
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The standard solution for tliis formulation using the least sqiiaies method 
(Bomford, 1980) gives: 

x„ = {a'^\va)~^a'^'wb 

V, = AX, -B 

Where the subscript a indicates that the value is approximate. The best sol- 
ution is obtained by iterating until X converges to less than some limit. 

W is the weight matrix. The standard daviation a, of a range measurement 
is an expresion of its accuracy. It is standard to weight the range measurements 
using which ensures achieving the minimum variance solution [Hamilton. 
1964], 

Assuming that there is no correlation between measurements, the weight 
matrix becomes diagonal as follows: 



l/aj 

0 

0 



0 

\la\ 



0 



0 

0 




Motorola estimates the standard deviation (1 o ) of a single range measure- 
ment to be 2 meters. The individual user has to set his own estimated ( o ) for 
his equipment (see section IV-H). 

The unbiased estimate of the variance-covariance matrix of the adjusted po- 
sition is given by 

Z, = oI{a''\VA)-' 



This matrix includes the variances and the covariances in the coordinate sys- 
tem we use ( X,Y coordinates) and appears as follows: 
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V _ 



(18) 




xy 



Assuming that a, = 2ni, then the combination of two ranges can result in po- 
sitional error of 7.8 m (drms) in the limiting case of 30® intersection angle. Drms 
is derived from the unbiased estimate of the variance covariance matrix as fol- 
lows: 

\drnis = ^'^((7^ -f Gy) 

This error of 7.8 m can be larger when the system is not calibrated. For un- 
calibrated systems the error of a single measurement can exceed 10 m [Motorola 
Mini-Ranger manual]. 

The quantity a; is the a posteriori variance of unit weight. It can be computed 
as follows: 



<iof 

Where dof is the degree of freedom, or n-2 when solving for X.Y position co- 
ordinates using n ranges. Ideally, crj should equal unity. 

According to Cross [1981]. when the degrees of freedom are small it is dan- 
gerous to use a single determination of a„ to judge the quality of the position. It 
is better to use an average \ alue of g^ from the whole data set. 

But the big question is "how we can say that the value of g^ is out of bounds? 
For this question a chi square one or two tailed test is the best indication, even 
though there are situations where it is not sufficient [Uotila, 1975]. If the test 
fails, this may be an indication that one or more of the following error sources 
exist that cause the large t’alue of g„. 

• Error in the mathematical model intolved in the solution (Scale factors. 

Refracthe Index used etc.) 

• Computational errors. 
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• High correlations between observations or v'ery poor geometry of the situ- 
ation. 

• Influence of the omitted higher order terms in Taylor series expansion which 
can be critical for higher order accuracy systems (primarily when the Master 
station is far away from to the Reference station). 

• Incorrect variance-covariance matrix of the observed quantities. It means 
that for a non well calibrated system, large values of are expected. 

• Blunders in observations arising from various factors in the measurement 
process (Multipath, unstable readings, etc). 

So if the statistical test used for cr^ shows that something is wrong an investi- 
gation is needed for the detection of the possible error source and its elimination 
from the data. 

B. PROPAGATION EFFECTS 

The line of sight limitation is one of the critical factors in the network instal- 
lation. The radius of curvature of the signal path is about 4 times greater than the 
Earth's curvature, so the heights of the Reference stations and the height of the 
Master station antenna govern the maximum range of the system, according to 
the formula. [Motorola Mini- Ranger Manual. 1981]: 

d — 4.04 (^^ /?,„ ) 

Where /?, and h,„ are the heights in the Reference and Master stations respectively 
and d the maximum range in kilometers. 

The maximum range of the system can be reduced due to the atmospheric 
conditions. Attenuation of the signal due to rain is given by [Casey. 1982]: 

A — ar^ 

Where A is the attenuation of the signal in db, r is the rainfall rate in mm hr and 
a.b are frequency functions. For the Miniranger mean frequency (5.5 GHz), their 
numerical \ alties are a= 1.48 x 10"^ and b= 1.1469. 

Another formula is gi\en for fog. based on the water concentration per cubic 



meter: 



A = 4.87a- 10"^V//Vz>//:w. 



Where M is the water content in grjm^ and f is the frequency in Ghz. 

Table 2 gives the attenuation of the signal for rain and fog. at 5.5 Ghz which 
is an average operating frequency for the Mini-ranger. In terms of distance, we 
have 8.0 km range reduction at 40 km for each db drop. So, for example, in a 
distance of 40 km with hea\y rain conditions we have 1 .2 db drop in the signal 
(1 1.6 km reduction in maximum range). This drop in signal power can introduce 
unexpected range errors when signal strength drops below the Critical Strength 
Threshold (CST). 



Table 2. ATTENUATION OF THE SIGNAL AT 5.5 MHZ 



Precipitation 
Rate (mm hr) 


Attenuation 

(db/km). 


Fog (gr m-) 


Attenuation 
(db km) 


2 


0.0015 


1.0 


0.0147 


4 (drizzle) 


0.0072 


1.2 


0.0177 


6 


0.0116 


1.4 


0.0206 


8 


0.0161 


1.6 


0.0236 


10 


0.0208 


1.8 


0.0265 


12 


0.0256 


2.0 


0.0295 


14 


0.0305 


2.2 


0.0324 


16 (hea\-}') 


0.0356 


2.4 


0.0354 


18 


0.0407 


2.6 


0.0383 


20 


0.0460 


2.8 


0.0413 



A possible problem that can appear in the sur\ey area is the multipath phe- 
nomenon of the signal. It can cause constructive or destructi\e inteference at the 
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receiver. There exist two types of multipath. The first, vertical multipath is 
caused by reflection of the signal on the sea surface. The destructive interference 
at the receiver due to vertical multipath is called "range hole". The second one 
is the horizontal multipath due to reflection of the signal on walls or other sur- 
faces. The destructive interference of the horizontal multipath is seen as "unstable 
readings". 

Range holes are very important in networks that operate at line of sight and 
can be a major source of positional problems. The best way to avoid such a 
problem is to put an additional Mini-Ranger Master station at a different height 
on the vessel. 

The accuracy of the range measurements is affected by the weather conditions 
(temperature, humidity, pressure). A fi.xed value of the refracth e index (N) is in- 
put into the system prior to a survey and kept constant throughout the sur\ey. 
The default refractive index (N = 320) corresponds to atmospheric conditions of 
20^C. 30 '’ 0 and 1013 mb respectively. 

Extreme atmospheric conditions, like high temperature along with either high 
relati\e humidity or dry atmosphere, can give rise to large scale errors in the 
measured distances. For this reason further investigation was done to estimate 
the magnitute of the atmospherically induced errors (see section IV-E). 

C. LOW SIGNAL STRENGTH 

Se\eral intestigators [e.g Casey, 1982] ha\e shown that range errors, both 
bias and standard deviation, tend to increase once the recei\ed signal strength 
drops below some critical level. This le\el, usually refered to as the Critical Signal 
Threshold (CST). is defined by the users accuracy requirements. Typical choices 
are between signal strengths of 8 to 18. The Motorola Falcon equipment itself 
flags signal strengths less than or equal to 13 as potentially suspect. 

The signal strength from the shore station is also a factor that can introduce 
large errors in the range measurements, but it is always output by the system and 
can be used in subsequent data processing. When signal strength is abo\e the 
CST limit, the errors ha\e a Gaussian distribution, (mean \alue zero and rela- 
tively small \ariance). But when the signal strength is below the CST limit the 
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Signal Strength 

Figure 3. VARIATION OF MEAN ERROR WITH SIGNAL SI RENGTM 

errors follow a logarithmic curve. That means the range error increases logarith- 
mically as one recei\ es smaller values of signal strength [Casey. 1982]. 

Figure 3 deri\ed by Casey, gives the CST at SS = 8 . For particular applica- 
tions each user can set his own estimated N'alue for CST in order to either reject 
the low signal strength meausurements, or to increase the standard dexialion of 
the measurement. 

The SS of each signal depends strongly on the distance from the reference 
station, (it drops with the imerse square of the distance). It can also wary due not 
only to multipath phenomenon, (as we explained earlier), but also due to the an- 
gle to the main lobe direction and the atmospheric conditions. 
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IV. TESTS ON THE MINI-RANGER NETWORK IN MONTEREY BAY 



A. MINI-RANGER CRUISE DATA 

On 31 March 1989 the Vessel Point Lobos, from MBARI, sailed for a trip in 
Monterey Bay to make trials and calibration on the Mini-Ranger stations. The 
Master station was installed on the mast at height 10 m above the water. During 
this trip useful data were collected and position estimates calculated. This data 
was processed in a preliminary fashion in order to detect possible error sources. 

Figure 4 shows the routes of the vessel followed. A div ision of the Bay into 
smaller areas was used, (as shown on the figure), in order to subdivide the data 
for processing and error analysis purposes. The division was based on the con- 
cept of nearshore and open b.';y regions. 

Another trip was subsequently done for NPS research purposes on 22 Sep- 
tember 1989, and some more data collected as shown on the same figure. The 
master station for this cruise was different and the data w collected under differ- 
ent conditions of atmosphere, time and instrumentation. It was considered possi- 
ble that the positional results from these two data sets may vary in accuracy. 

First we e.xamined the earlier data set of MBARI. which already had the 
positions and their statistics calculated. We were informed by MB.ARl that for 
the position derivation they used the method of \ ariation of coordinates. No de- 
tails on the specific code in this software were avaliable. 

The first step was a general view on the residuals and the "standard devi- 
ation" provided in M BAR Is data of the whole trip. This showed that something 
did not work vv^ell when the vessel was nearshore , especially in the SW portion 
of the Bay. Figure 5 shows the areas of the greater and more v'ariable standard 
deviation of the data set. In general as shown in Figure 5. areas E and G have 
fewer problems. It was not obvious which station or stations combination was 
responsible for the big residuals. 

.Another problem was the v'ariation of signal strength and its unexpectedly 
high values. In many portions of the data set, the signal strength was much higher 
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Figure 4. SHIP ROUTES FOR IFIE TWO DATA SEIS/AREA DIMSION 
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Figure 5. AREAS OF GREATER STANDARD DEVIATION 
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than expected for such a range. The appearance of signal strengths greater than 
70 in distances more than 10000 meters is unusual. This is a strong indication of 
possible system malfunction. 

Because the data set did not show the expected accuracy, the following series 
of tests and checks were performed: 

• A visit to the stations for inspection of directional limitations, antenna in- 
stallation problems and other possible affecting factors. 

• A check of the station positions with geodetic survey techniques, (triangu- 
lation, traverse, or differential GPS), and a least squares adjustment of the 
survey data with fixed geodetic control points around the Bay. 

• Check of the software used for the position derivation using the variation of 
coordinates algorithm with accurate grid scale and slope distance correction 
techniques. 

• A theoretical view of the geometrical accuracies with the use of various sta- 
tion combinations in the Bay area. 

• A check on the effects of extreme meteorological conditions to check for the 
possibility of introducing range errors in e.xcess of 1 meter. 

• A range hole prediction for the network. 

• An investigation into the possibility of correcting those ranges with signal 
strengths below the CST. 

• A calibration test for determining biased or linearly varying errors and their 
\‘ariability with signal strength abo\ e the CST. 

B. INSTALLATION CHECK AND ANTENNA DIRECTIVITIES 

On 20th of July, a visit to the stations TREVOR. WATS and PACKARD 
was done for \isual examination, check on the antenna directi\ity , possible sig- 
nal blocking from nearby obstructions and other possible affecting factors. The 
stations HA'i'S and HANK were \ isited on the 1st of August and 25th of Sep- 
tember respectively. The alignments of the antennas were determined to approx- 
imately 2 ° using a small portable magnetic compass. All the magnetic bearings 
are given as true bearings. A short description is gi\en below for each station and 
a map is presented (Figure 6) for clarification; 

• TRE\'OR: This station is installed on the roof of a house in Santa Cruz. The 
antenna is attached to the chimney of the house and is oriented towards 
175^. The signal is blocked by a rocky hill from 215® to approximately 100® 
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Figure 6. OBSERVED ANTENNA DIRECTIVITIES OF MINI-RANGER 
NETWORK IN MONTEREY BAY 
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(clockwise). There is no clear view of the Northern near shore Bay region but 
this area is of little interest. 

• PACKARD: Installed on the roof of a wood craft building, hooked on a thin 
mast, on private property between Watson\’ille and Moss Landing. Its an- 
tenna points to 250'’ and may be blocked by the Moss Landing electric plant 
stacks between 202'’ and 205" There were no other problems observed for this 
station. 

• WATS: Installed on the roof of a commercial building near the beach at 
Watsonville. Its antenna points at 215" .The entire Bay area is unobstructed. 
The antenna is approximately 2 feet from the \’ertical wall of the building. 

• HAYS: Installed on the roof of Silas B. Hays hospital at Fort Ord with an 
uninterupted Bay view and directivity of its antenna to 295" . The perfect 
\ iew along with the height of the station makes it useful in some portions 
SW of Point Pinos . The signal is blocked from the higher Pacific Grove hills 
from 253" back to the South. 

• HANK: Installed on the roof of a house in Pacific Grove. While it has good 
height the signal cannot reach the nearshore portion from Point Pinos to 
Monterey harbour due to high trees existing downhill. The directh’iiy of its 
antenna is 345" and the signal may been blocked by high trees from 005" to 
035". The station is offset from the position specified by the preliminary co- 
ordinates. This particular correction is given in the discussion in section C 
below. 

• DOPPLER: This station was installed near NPS Beach Lab. but due to 
safety reasons was subsequenly removed. 

It was not possible to check the electric current supply of the stations, but a 
\oltage stabilizer can help the system to be kept healthier and unstable outputs 
can be a\oided in the future. 

C. POSITION CHECK 

Geographic coordinates for the station sites had been determined on the 
North American Datum of 1983 (NAD S3) by geodetic sur\'ey techniques, with 
ties to existing horizontal control points published by the National Geodetic Sur- 
vey [Schnebele, 1989]. Conventional ground sur\ey and Global Positioning Sys- 
tem (GPS) techniques were both used. The sur\'ey data had been analyzed with 
a three-dimentional least-squares adjustment software package (GEOLAB, by 
GEOSurw Inc.). Absolute horizontal position accuracies estimated by the ad- 
justment were better than 8 cm (circular error. 95% confidence level). Additional 
work to confirm the offset positions of the Mini-Ranger units at P.ACKARD and 



WATS, has yci to be completed, but is not expected to shift their preliminary 
positions by more than a few centimeters at most. 

These station coordinates were con\crted to Universal Transverse Mercator 
grid coordinates (UTM Zone 10) for use in the vessel position computation algo- 
rithms. The UTM system has been adopted as the standard grid used by MBARl 
and NFS with this network. The specific Transverse Mereator Projection algo- 
rithm [Floyd, 1985] is accurate to within 1 cm for the entire zone. 

As described earlier, station HANK was offset from the preliminary position. 
The offset was measured by a simple magnetic compass for azimuth and a tape 
for distance and should have errors no greater than 5 cm. The measured bearing 
was 066® true and the distance 41 inches. The two positions were at the same 
height, so the corrections in UTM arc: dx = 0.952, dy = 0.424, dz = 0.00 in meters. 

The results are sumarized and the revised station positions provided in the 
Table 3. With the e.xeeption of station HANK, the revised positions of the 
stations differ by less than 10 cm from the preliminary positions. It is concluded 
that the station positions are not the source of the obser\ ed inaccuracies. 

D. ESTIMATION AND COMPARISON OF DERIVED POSITIONS 
1. Development of position derivation algorithm 

The measured ranges were used to recompute the ship's positions by the 
method of \ ariation of coordinates [Cross, 1981]. These recomputed positions 
were then compared with those provided by MBARl. The comparisons showed 
that there is a discrepancy between the given positions and the recomputed posi- 
tions on the first data set. which in some cases were significant. Because of these 
discrepancies, we ignored all the git en MBARl positions and recomputed all their 
data. 

For this purpose, a program was written in Turbo Basie. This program 
uses the \ ariation of coordinates method on the UTM projection surface. Scale 
factors were computed for each line. The program is given in the Appendix 1 and 
can be followed by comments pro\'idcd in the code. 

The measured distances from the system, for each case, are corrected for 
their slope with the assumption of a right triangle (ignoring the eur\ ature of the 



Table 3. REVISED POSITIONS OF MINI-RANGER 



STATIONS 



Name 


Latitude(N) 
Longitude(W) 
deg min sec 


UTM coordinates 


Elev^ation 


Trevor 


36 58 28.785 
122 02 31.536 


585260.158 

4092490.287 


51.1 


Wats 


36 54 26.896 
121 50 39.821 


602947.636 

4085231.429 


22.6 


Packard 


36 49 44.404 
121 46 04.911 


609863.147 

4076611.241 


33.4 


Hays 


36 38 33.824 
121 47 45.894 


607621.298 

4055915.205 


136.7 


Hank 


36 36 25.393 
121 55 08.683 


596670.450 

4051826.796 


134.1 



Earth). This assumption inserts an error of 4 to 5 cm at the extreme ranges which 
is insignificant for Mini-ranger measurements in the areas eovered by the net- 
work. but it makes the calculations simpler. For the complete correction one can 
look in chapter 15 of Laurila [1981]. 

The next correction that was made to the measured distances was 
for the scale factor of the UTM projection. The following formula computes the 
scale factor for the line between the Master and each Reference station: 



Where m is the required scale factor, is the scale at the central meridian 
(0.9996), R is the mean radius of the Earth (6375 km) and is the mean Easting 
N'alue for each case. The mean Easting ^'alue is derived from the formula; 

El = + £,£, + El 

Where £, and E 2 are the true Eastings (x, — 500000), of the Master and the Ref- 
erence stations respectively. The formula is a simplification of the one given by 
Bomford [1980, pi 93] and the error due to this is of millimeter order, which for 
our purposes is negligible. 

The outputs of the program are the least squares estimate for the position 
of the \essel, the variance co\ ariance matrix the residual \ector (V) and the 
a posteriori t'ariance of unit weight al. Measurements were assumed independent 
so the diagonal weight matrix is treated as a vector in the calculations. 

The a posteriori variance of unit weight al is the one that is used for as- 
sessing both the accuracy of the network and the ranging data. ol can be scaled 
as is well known, by increasing the standard de\iation of a single range meas- 
urement. The resulting n; is examined statistically with the chi squared, two 
tailed test, at a 95*'/o confidence le\el. The null hypothesis. ol=\, was tested 
against the alternatite When the null hypothesis is accepted, it is an in- 

dication. assuming no data blunders, that the system works within the specifica- 
tions pro\ ided by the manufacturers. When the null hypothesis is rejected se\ eral 
possible reasons for this exist. These were summarized in section Ill-A. 

It also should be noted. howe\er, that the presence of any unresohed 
systematic effects tends to inv alidate this tv'pe of test. One potential source, not 
modelled in the position computation algorithm, involves the non-simultaneity of 
range measurements on a moving ship. In effect, the antenna is not stationary 
over a cycle of range measurements to all Reference stations. Given each range 
measurement takes about 55 milliseconds, four range measurements may be 
spread over a period of about 0.15 seconds. At a nominal v'essel speed of 6 knots, 
the possible error is about 0.45 cm. which is of little importance in this anah’sis. 
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Additional error can be inserted due to mast motion. In rough sea condi- 
tions the error is big enough (more than 1 meter) and it has to be taken into ac- 
count by increasing the estimated standard deviation of the range measurements. 

2. First data set 

The results of the chi square test at a 95% confidence level using standard 
de\'iations of 2 and 3 meters on the measured ranges respectix’ely, appear in Table 

4. 

The acceptable values of cl under these conditions are: 

For three ranges 0.00 < cl < 5.024 

For four ranges 0.025 < cl < 3.69 

The numbers gi\ en in the first two columns of Table 4 represent the per- 
centage of the number of positions for which the null hypothesis was rejected 
against the total number of the positions examined. The number pro\ided does 
not include some abnormaly high a posteriori variances of unit weight which were 
caused by multipath effects, or low SS. These data were rejected as blunders. 
Their percentage appears in the last column. 

As can be seen. e\en when c =3m. the number of rejections is clearly 
much higher than the 5'’ o which one would normally expect. 

In the column 'Difference' we can see how many of the positions that 
were rejected when c =2m. arc accepted with c =3m. This is an indication ei- 
ther of poor calibration of the system or that the system does not perform to the 
standard claimed by the manufacturer. It should be noted that the range data 
had been corrected for an estimate of systematic errors determined in a prior 
calibration. 

The results shown in Table 4 suggest the possibility that the master sta- 
tion used may ha\e beeen poorly calibrated. For this reason, an examination of 
the distribution of the residuals for each reference station was undertaken. After 
examining the distribution of the range residuals for each reference station, it 
became clear that some stations were consistently biasing the position solutions. 
From the distribution of these residuals, the bias was estimated by comparing 
their mean, median and most probable values. The mean residuals and the 
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Table 4. REJECTED PERCENTAGE WITHOUT REVISED 
CORRECTORS 



AREA 


0 = 2 


0 = 3 


Difference 


Rejected 

blunders 


A 


68% 


48% 


20% 


11% 


B 


82% 


37% 


45% 


9% 


E 


69% 


22% 


47% 


25% 


F 


27% 


13% 


14 % 


33 "'o 


G 


43<’'o 


27% 


16T'o 


29% 



adopted correctors are shown in Table 5. This simple averaging of obsert ed res- 
iduals is not an ideal procedure for estimating bias errors. The least-squares 
process tends to redistribute an error in any one range to errors (residuals) in all 
ranges. Recognizing this limitation, however, the procedure gives a useful esti- 
mate of uncorrected bias errors. 

The results for the corrected data are shown in Table 6 and can be di- 
rectly compared with Table 4. The improvement as we see is above 65 To. 

In the column 'Difference' we can see. in comparison with the one of table 
4. that the \ alues are much lower, which is an indication that most of the cali- 
bration uncertainty has been removed. 

Even after this procedure, some a posteriori variances of unit weight are 
still high, especially in areas A and B. The possible reasons for this can be: 

• Incomplete removal of the blunders. 

• Operation of some stations at low SS which in turn causes larger deviations 
from the normal values. 

• Incomplete removal of the bias errors described above. The technique used 
does not eliminate the entire error, but only estimates it. 

• For any one of a number of reasons the signal can be lost temporarily. 
When it is regained it sometimes lacks stability for several readings. 
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Table 5. MEAN RESIDUAL AND VALUES USED FOR 



RANGE CORRECTORS 



STATIONS 


Mean residual 


Used corrector 


Trevor 


3.91 m 


-3.5 m 


Wats 


-1.66 m 


1.5 m 


Packard 


-0.24 m 


0.0 m 


Hays 


2.61 m 


-2.5 m 


Doppler 


3.64 m 


-3.5 m 


Hank 


-0.24 m 


0.0 m 



Table 6. REJECTED PERCENTAGE AND IMPROVEMENT WITH RE- 



VISED CORRECTORS 



AREA 


( 7=2 


(7 = 3 


Difference 


ImproN'e- 

ment 


A 


38'*^o 


227 « 


16% 


44'* 0 


B 


31% 


14«o 


17% 


62 % 


E 


14^0 


06G) 




79 o 


F 


8''o 


O'^o 


8% 


70'^ 0 


G 


13‘’o 


l^/o 


6**0 


69®/o 



• Higher \ alue of the standard deviation of the system used than the one used 

for the weight matrix formation. 
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3. Second data set 



The second data set was examined statistically in the same way and the 
results appear in Table 7. We can see that the second data set gives much better 
solutions. This was expected because it was known that the master station used 
for this data collection, was well calibrated. 



Table 7. REJECTED PERCENTAGE FOR THE SECOND 
DATA SET 



AREA 


Rejected percentage for c = 2 


A 


9% 


B 


6^0 


C 


6*'/o 


D 


7% 


E 


7^0 


F 


7% 



The previous data set. for reasons noted earlier, showed major problems 
when the ship was in areas A and B. The second data set. howet er, shows con- 
sistent results in all areas and a low rejection percentage even when o =2m. 
These results suggest that the standard deviation of the range measurement 
should be between 2 and 3 meters. 

E. GEOMETRICAL ACCURACIES 

The position quality given from the Mini-ranger network in Monterey Bay 
depends strongly on the relatit e geometry of the ship and the reference stations. 
The geometric position quality can be represented in terms of error ellipses, con- 
fidence circles, or with several other appropriate techniques. When several 
stations are used, computer techniques are necessary if detailed analysis is to be 
undertaken. Such techniques are described in Perugini [1988]. 
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The provided figures in Appendix 2 were taken from computer calculations 
using software associated with MBARI's Navigation processor. The numbers 
pro\'ided represent the semi-major axis of the predicted error ellipse at 95% con- 
fidence level that is formed at the particular position, assuming the standard de- 
viation of a single range measurement is 3 meters. 

The best working combination appears in Figure 7 and is obtained using 
stations TREVOR, PACKARD, HAYES and HANK. We can see that the ge- 
ometry of the position is very good with two combinations of only three stations, 
i.e; station combinations TREVOR, PACKARD, HANK and TREVOR, 
P.ACKARD, HAYES (see Appendix 2) 

As is easily seen, station TREVOR occupies a position that gi\es strong ge- 
ometry in the mid-Bay regions and its position is the most critical for the network. 
On the North side of Monterey Bay there are also other higher positions on which 
a permanent station could be installed with a much better view of the Bay. Most 
of them are in the Uni\ersity of California Santa Cruz campus and some 
arrangments ha\ e to be made. 

In the Watson\ille area the hills o\er Buena Vista Point have a good Bay 
\iew but the problem of a permanent electrical supply and security has to be 
sohed. Such a position would overcome the height problem of the WATS station 
while maintaining its geometry. 

Stations PACKARD and HAYS are at the highest ele\ations in their areas. 
There is no ad\antage to shifting their positions without deteriorating the geom- 
etry. 

Station DOPPLER does not add much to the geometric quality of the net- 
work in the mid-Bay but can be used for positioning near Monterey harbor where 
signals from PACK.^RD and HANK are obstructed. It also could be shifted NW 
of Point Pinos on the lighthouse or on another high position to impro\e the ge- 
ometry in the Monterey canyon area. Further in\estigation on this station has 
to be done for the best choice of a permanent position. 
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riguie 7. BEST WORKING COMBINATION OF THE MINI-RANGER 
STATIONS 
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Station HANK gi\ es very good geometry but its position has some limitations 
because of obstructions to its Bay view. Some station at the southern site of the 
Bay is necessary to maintain a good network geometry for the areas of interest. 

Generally we can say that the geometry of the network is not responsible for 
the big standard deviations observed in most regions around the Bay, but it be- 
comes more important as one moves offshore near the Northen and Southern re- 
gions. 

F. ATMOSPHERIC CONDITION EFFECT OF THE AREA ON THE 
MICROWAVE PROPAGATION VELOCITY 

An in\estigation was carried out to see the maximum possible error due to 
propagation \elocity which could occur under e.xtremc atmospheric Bay condi- 
tions as compared with mean Bay conditions. The atmospheric formulation is 
described in Chapter 6 p.l33 of Laurila [1981]. The effects of temperature (T). 
pressure (P) and humidity (E) on the refracti\e index (N) are given from the 
Essen formula: 



P n 9'> 10 

A = 77.62 Y -( - 27.19 -^ )£ 

Where P and E are in mb and T is the dry bulb temperature in K°. But 



£ = 



100 



Where RH is the relative humidity (%) and £^ is the saturation \apor pressure. 
gi\ en in mb. at the dry bulb temperature. £^ is a logarithmic function of temper- 
ature so the effeet on the refraeti\e index depends on the temperature and the 
relathe humidity as shown in Figure 8, for various relative humidities and eon- 
stant atmospherie pressure at 1013 mb. 

Daily weather logs, kept by the Meteorology Department of the NPS. were 
searched for the maxima and minima of temperature and relatixe humidity in 
Monterey during 1988-89. This was done so as to eompute the a\ erage x’ariability 
of the refractit e index and the possible error in distance that is introduced foi' the 
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Pressure= 1013 mb 

Hum.30% Hum.50% HLm70% Hum.90% 




Temperature (C) 



Figure 8. VARIATION OF THE REFRACTIVE INDEX 

actual atmospheric conditions. The base for the comparison is the ieliacii\e index 
\ alue (N = 320) that is used from the Mini-ranger for the range calculation. 

In Monterey Bay during the spring, temperature averages about 55 °F and 
humidity \ aries from SO’f'o during the hotest days to 95% duiing the cold nights. 
During working hours, a relatise humidity of 55% with tcmperatuies 70 P are 
usual. Extreme values of 95% RH and 55 <>F are rare, but if they are accompa- 
nied with high atmospheric pressure, 1030 mb, then they gi\e N = 3d3.6 which 
in turn corresponds to a distance error of 23.6 ppm. In a po.ssible distance of 
40000 m the introduced error is about I meter. On the other hand, low atmo- 
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spheric pressure of 1000 mb with 80 °F and 20% humidity, gives N = 288 and 
32 ppm error, which in distance terms means 1.28 m at 40000 m measured dis- 
tance. Average \ alues of N for the spring during working hours are between 310 
and 330, which indicates that the assumed Refractive Index of 320 is fair. 

During the summer the increased temperature is accompanied with a drop in 
the relative humidity so the result tends to be the same as in spring with an av- 
erage N = 320 to 325. The extreme cases here also rarely introduce errors more 
than 1 meter at distances of 40000 m. 

The fall season is the one with a slightly different mean N and more vari- 
ability. Warm weather with high humidity is more likely to occur. In the extreme 
cases, at maximum distance, the error could reach 2 meters, but it rarely occurs. 
The average N = 330 is better to use. 

Winter shifts the refractive Index back to the lower values due to the lower 
temperatures of the atmosphere, which makes the extreme cases more rare. A 
mean \’alue of N= 315 is better for use. 

Generally the resulting error should be less than a meter at extreme ranges 
and sev eral decimeters on average, when extreme atmospheric conditions happen. 
Even though the error is small, the use of the seasonal refractive index is advised. 

G. RANGE HOLES 

Range holes occur due to the vertical multipath, which is caused from re- 
flection of the signal from the sea surface. This is common to all microwav e sys- 
tems that use a line-of-sight propagation path. Figure 9 [Gilb-Weedon. 1976] 
illustrates the situation. 

We count the range holes from the first to the Nth. Here N represents the 
integer number of wavelengths difference between the direct and the reflected 
paths. The first range hole occurs at the longer distance from the shore station. 

Several factors affect the range hole characteristics: 

• The distance between the reference and the master stations. The longer the 
distance, the wider the range hole zone. 

• The height of the stations. The higher the reference stations the fewer the 
range holes near them. 
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• The wavelength of the signal. This affects both the width of the range holes, 
and the distance from the reference station at which the}' occur. 

• The atmospheric conditions along the signal path cause loss of SS. I hus the 
reflected signal Horn the sea surface, due to additional losses from the re- 
flection. is too weak to interfere the direct one, when they both arri\e to the 
master station. 

• The sea surface roughness. This affects the type of reflection which occurs. 
In a strongh' wa\e modulated sea, the appearance of the range holes is rare. 

• Changes in the sea letel due to tides, winds etc., cause range holes to shift 
location horizontally. 

Small reflection angles cause more frequent range holes because the signal has 
lesser losses in comparison with the case when the reflection angle is bigger. The 
additional losses of SS. in the second case, weaken the rellected signal vhich can 
not then interfere with the direct one. 

In reality, the destructive interference is rarely total due to losses of the re- 
flected signal from the sea surface, but it can reduce the SS to below the sensi- 
ti\ ity of the receiver or below the CST in the case of Mini-ranger. This will result 
in signal loss, or worse, in a bigger standard deviation of the measurement. 
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In the Monterey Bay network, the range hole predictions arc shown in Table 
8 for each station. For the maximum range of each station we assume a ship an- 
tenna of 10 m. 

Range holes closer than the third occur rarely because of the bigger angle of 
incidence which attenuates more of the reflected signal. 



Table 8. RANGE HOLE PREDICTION FOR MONTEREY BAY 



STATIONS 


Maxi- 

mum 

range 


Range for the 
first hole 


Range for the 
second hole 


Range for the 
third hole 


T revor 


42000 m 


18972 m 


9486 m 


6324 m 


Wats 


30000 m 


8477 m 


4238 m 


2826 m 


Packard 


37000 m 


12440 m 


6220 m 


4147 m 


Hays 


61000 m 


50348 m 


25174 m 


16783 m 


Doppler 


26000 m 


3670 m 


1835 m 


1223 m 


Hank 


60000 m 


49358 m 


24679 m 


16453 m 



Table 8 is based on the following formula [Gilb-Weedon. 1976]. in which a 
plane earth is assumed. R represents the distance from the reference station to 
the range hole center in meters. 



Ih 

^ = -">77 



The distance of the maximum expected signal is: 



/h 



= 4/h. 

‘ {2k + 1) 



Where / 2 , and lu are the heights in meters of the reference and the master stations 
respecti\cly. 1 is the wa\clength and k is an integer (k = 1.2.3....). 
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1st Hole 2nd Hole 3id Hole 




Figure 10. RANGE HOLE SHIFT VS 1 IDAL HEIGHT 

In the Bay iheic is a sea level variation due to tidal effects of 6 to 7 feet 
maximum. This effect tends to shift the range hole center towards the shore sta- 
tion at high water and away from the shore station at low water. The o\erall ef- 
fect of this phenomenon for the stations appear in Figure 10, for each range hole. 

Typical wind and wave conditions in the Bay are such that 'perfect' signal 
reflections are not expected. Thus range hole problems should be rare and limited 
to conditions of calm sea accompanied by light winds. Range hole plienomena 
were not detected in the positioning data studied. 
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H. CRITICAL SIGNAL THRESHOLD (CST) DERIVATION 

The relatively large distances and varying propagation conditions encount- 
ered in the Monterey Bay network result in many low signal strength readings. 
It is important that a method be found for using these readings if at all possible. 
Casey, 1982. has shown that low signal strength ranges can be used if corrected 
in postprocessing by an empirically derived relation between bias error, standard 
deviation, and signal strength. 

The recommended procedure to establish this relation is to perform repeated 
range measurements over a known length baseline, with a variable attenuator to 
artificially vary the recei\’ed signal strength. These attenuators were not 
at aliable. but the effect was achie\ ed by ranging over relatively long baselines of 
several different lengths, in order to obtain readings over a spread of low and 
high signal strengths. Two sets of range readings were collected for different 
combinations of Master and reference stations, with signal strengths varying from 
about 5 to 25. 

For each Master Reference station pair, ranges with signal strength o\'er 15 
were a\eraged to estimate a mean range, free of low SS effects. Differences be- 
tween this mean range and all low SS ( < 15) ranges were computed, representing 
the additional error due to low SS. By this technique, error estimates from dif- 
ferent Master Reference station pairs could be combined without concern about 
systematic differences between equipment pairs. The resulting range error esti- 
mates were grouped by signal strength; 5-6, 7-8, 9-10 and so forth. Mean errors 
and standard detiations were then computed for each signal strength group as 
shown in Table 9. As expected, both the mean error and standard de\ iation tend 
to increase as the signal strength decreases. 

The two data sets were compared to ensure repeatability of the results before 
deriving the empirical relation. In Table 9 we can see the results of the statistical 
tests for the two data sets. They are treated as samples from one population, but 
each with their own estimated means and standard det iations and examined with 
the Fisher - Behrens statistical test [Hamilton. 1964] at a 95*’o confidence le\el. 
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as follows. The T test limit was taken from the Student's t distribution table with 
degrees of freedom calculated by: 



dof= 



2 2 
(7 1 ^ 1 ') 

Hi ^ !U ’ 
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+ 
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02 



U2>2- U 



The calculated value for the test is gi\en by: 

(.Y, - J2) - (/M - t-h) 



Where aj and 0 ^ are the estimated variances, //, and jUj are the estimated mean 
errors, and /i, and //j are the number samples in the first and the second data sets 
respectively. 

We examine the hypothesis: 

//o :.Y| — .Y 2 = 0 The two data sets belong to the same population. 

Against the alternati\e 

Hi'.Xi — X 2 ^ 0 The two data sets don't belong to the same population. Table 
9 shows the deri\ed statistics from the two data sets and the limit of the T test in 
each case. E\ery calculatetl value, for the Fisher-Behrens statistical test, that is 
bigger than the T \ alue means that the hypothesis is rejected so something else 
going on in the data sets. 

In the column 'Calculated value', the number is deri\ed using the means and 
the standard deviation of the measurements. We can sec that when SS is 5-6 and 
17-18, the calculated value is bigger than the T test limit, which means that in 
these cases the null hypothesis is rejected. 

The number gi\en in the column 'Calculated \alue with increased 0 ' is a 
combination of the standard de\ iation of the long term variation in the instru- 
ment and its observed standard dexiation. As we can .see the result is a total 
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Table 9. STATISTICAL TESTS FOR THE I WO SAMPLE DATA SEIS 



s.s 


1st data set 
(mean error / a) 


2nd data set 
(mean error/ o ) 


T test 
limit 


Calcu- 
lated 
value (u) 


Calcu- 

lated 

value xvith 
increased 
o 


5-6 


15.1 4.5 


13.5,5.3 


1.98 


2.48 


1.78 


7-8 


8. 1/3.7 


7.8, 3.8 


1.98 


0.63 


0.41 


9-10 


4. 5/3.0 


3.8'2.6 


1.98 


1.68 


0.98 


1 1-12 


1.9 1.7 


1.9/ 1.7 


1.98 


0.28 


0.19 


13-14 


0.3 1.8 


0.3, 0.7 


1.98 


0.28 


0.2 


15-16 


-0.0 1 .5 


0. 1/0.6 


1.98 


0.98 


0.37 


17-18 


0.8/ 1.4 


0.3, 0.7 


1.98 


3.35 


1.17 



statistical agreement of the data sets. We conclude, therefore, that they belong to 
the same population. 

The ne.xt step is the derivation of the equation that will permit the use of the 
lower SS more effeetixely. The equation derived was 

( - 7 . I -P . 2 1 -27.206) 

Inx 

Where y represents the mean error in meters and x the low signal strength 
and the given numbers represent meters. The resulted y \ alue represents the ob- 
serxed minus the calculated (C-0) range correction. 

Figure 13 shoxvs this resulting best fitted curxe. The correlation for the case 
is more than 0.999 xvhich indicates that the resulting equation is xery close to 
perfectly modelling the data. This equation can be included in the software as a 
correction for the loxver SS rather than rejecting them bcloxv the CST. 
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Figure 1 1. FFFTED CURN E 70 TFIE COMBINED DATA FOR TESTS 1 AND 
2 

Additionally, due to the increased standard de\ iation of the ranges at lower 
SS a higher value for standard de\iation of a single measurement should be used 
in the weight matr.x formulation. 

1. V ARIATION OF THE STANDARD DEVIAT ION WITH TIME 

Eigth sets of range measurements were collected o\er the known distance be- 
tween station HAYS and geodetic station RANGE 7. The choosen distance was 
short (2341 m) to eliminate the effect on distance of the \ariation of weather 
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conditions. These 8 sets were collected on 8 different days and they have length 
of about 1 hour, with 3 second intervals between range measurements. Table 10 
gi\ es the bias of these measurements with respect to the known distance and the 
standard de\ iation of a single measurement. The results from Table 10 show that 
the system has a standard deviation which is indicative of its precision rather 
than its accuracy. In addition, it has a variable bias. The total standard deviation 
for the sy'stem must reflect both these components. 



Table 10. STANDARD DEVIATION AND 
BIAS OF RANGE ERROR 



NUM 


Standard dex i- 
ation 


Bias 


1 


0.285 


-0.725 


0 


0.206 


1.447 


3 


0.255 


1.324 


4 


0.266 


1.531 


5 


0.374 


0.577 


6 


0.856 


0.007 


7 


0.576 


0.501 


8 


0.541 


1.173 



The biases shown in Table 10 were used statistically to show that the stand- 
ard deviation of one measurement is 2.54 meters, which concurs with the \ alue 
deri\ ed by evaluation of aj estimates in section C of this chapter. The procedure 
followed for this derit ation was a simple multiplication of the standard dex iation 
of the mean bias found from the data set by the square root of the number of 
data sets. The user shouldbe aware that the standard deviation can be larger on 
board a xessel in a rough sea due to mast moxements and the selcctixity of the 
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stations that the system makes(non simultaneous measurements to the reference 
stations). The bias variation from day to day makes the calibration of the 
stations time consuming because a one day calibration will not remo\e all the 
bias. 



Table 11. MEAN DIFFERENCES AND STANDARD DEVIATIONS FOR 
STATIONS HAYS AND TREVOR 



DIS- 

TANCE 


Mean differ- 
ence (Hays) 


Standard devi- 
ation (Hays) 


Mean differ- 
ence (Trevor) 


Standard de- 
viation 
(Trevor) 


SHORT 


0.55 


0.8 


- 


- 


ME- 

DIUM 


1.02 


0.55 


0.91 


0.48 


LONG 


1.27 


0.75 


1.48 


0.68 


V.LONG 


1.2S 


1.2 


0.84 


0.95 



J. VARIATION OF THE STANDARD DEVIATION WITH DISTANCE 

The second thing was the variation of the standard deviation with the dis- 
tance or with S.S. when it is abo\e the CST. Se\eral problems limited the data 
collection, like the unavailability of an attenuator, the poor network of known 
positions around the Ba\' and the low ele\ ation of possible calibration sites. The 
data were sufficient only for the exaluation of station HAYS and partially of 
station TREVOR. 

Table 1 1 shows the mean difference of the measured distances from the true 
distances and the standard deviations of each set of range measurements. We can 
sec that there is variation of the mean error with distance in both stations but no 
indication e.xists that a function can be dcrixed because the mean differences with 
the gix en standard dex'iations are not statistically different. Actually, a functional 
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variation of the standard deviation with distance was not expected, when oper- 
ating above CST. 
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V. DISCUSSION AND CONCLUSIONS 



The Mini-ranger network in the Monterey Bay was installed for research 
purposes by MBARI and NFS and was examined from various aspects. Useful 
results were drawn that can help in the understanding of the operation of trans- 
ponder type systems in other areas, and that will make the present network more 
useful when operating at low SS. 

Figure 12 shows the o\ erall coverage of the Mini-ranger with all the limita- 
tions taken into acount. except that weather conditions can attenuate the signal 
and in turn reduce the maximum range of the stations. For the maximum range 
derivation a height of the master station on the ship was assumed to be 10 meters. 

From this network and the data processed, the following conclusions can be 
drawn; 

• For the calibration of such systems, if one wants the maximum accuracy that 
can be provided, a bias taken from single obsert'ations over known distances 
or with the baseline method may be not enough. It is better to use repeated 
obser\ ations and a mean corrector derit’ation. 

• The positions of the stations in the Bay were, with the exception of station 
HANK, determined to a sufficint let el of accuracy such that the\ were not 
responsible for some of the large range standard deviations. 

• The software of the user is important in the position derivation and can 
cause unexpected errors if not properly validated. 

• The meteorological conditions in the Bay can cause attenuation of the signal 
and. in turn, reduce the maximum range ability. The range error from the 
delay of the signal due to meteorological conditions does not exceed the 1.2 
meters in the extreme conditions and is smaller than 0.2 meters under the 
usual operational conditions. 

• Range holes in the Bay can occur, but should be rare due to usual windy 
Bay conditions and the modulated surface from the swell. 

• We ha\e shown that ranges measured at low SS need not be rejected, but 
can be corrected using an appropriate error function. This function may vary 
from system to s\'stem. 

• We show that the standard det iation can be considered to consist of two 
components, the unprcdicted bias, which is smaller when the ststem is well 
calibrated and the measured standard deviation of the data. Generally the 
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Figure 12. MINI-RANGER COVERAGE IN MONTEREY BAY 



a of one measurement under ideal conditions is a little greater than 2 meters. 
It is more proper to use a = 3 meters when the system operates at sea. 
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There are no indications for variation of the standard deviation with dis- 
tance (or SS) above CTL. 



APPENDIX A. COMPUTER PROGRAM LISTING 



REM PROGRAM "LESS" FILE; NICK3. BAS 

NICK KRIONERITIS Date: 8-8-89 

Modified 10-17-89 by KJS 

VrVr Vc Vc V'T V? Vc Vr Vc Vc Vr Vc V? Vr Vc Vc 

REM This program calculates the optimum position of 
a ship using distances from known stations. It 
serves the solution for a miniranger network with 
accuracy better than 10 centimeter. 

REM The subrut ine SOLVERANGX was changed from CDR K. SCHNEBELE 
program "HYDROPLOT". 

V* V^V^^rVr Vc Vc Vc Vc Vc%Wc Vr V? VcVc V% Vc Vc Vc V'T Vc Vr Vc Vr Vc Vc Vc Vc Vc 



$STACK 9216 

$INCLUDE '^SOLRANGX. SUB” 

$INCLUDE ”NICKMR. SUB” 

ON ERROR GOTO ERRORHANDLER 

Print ”Enter the number of stations” 

Input ”N?o=”,N% 

Lprint ’’Number of stat ions=” ,N?o 



DECLARATION OF THE VARIABLES 

AP//(1): Approximate position of theship in meters (UTM) 
SP//(2): Station positions in 2D array in meters (UTM) 

R//(l): Ranges observed on the receiver in meters 
LX//(1): Stations Eastings for the scale factor affection 
X?'/(l)’ False Eastings of the stations 
Y 7 KI): Northings of the stations 
Z/r(l): Elevation s of the stations from MSL 
SF?'/(1): Scale factor calculated for each case 
W!(l): Diagonal Weight matrix assuming indepedence 
among the measured ranges 

C0DE%(1): Index numbers of the 4 stations in each record 
AllSta#(2): UTM Coords (x,y,z) of 6 stations possible to use 



Vectors for S0LVERANG2 but not used in this version 
Cm0s!(l): Bias corrections for each of 6 stations 
Vgts!(l): Weight factor for each of 6 stations 



DIM AP#(1: 3) 

DIM DYNAMIC SP#( 1: N%, 1: 3) , R:Kl:N%) , LX#(1:N%) 

DIM DYNAMIC X#(l:N?o) , Y#(1:N%) , Z//( 1: N%) 

DIM DYNAMIC M//( 1 :N?;) , SF//(1:N?0 , BRes#(l:N?0 , W!( 1 :N%) 

DIM code%(l: 4) , AllSta//( 1: 6 , 1: 3) ,CmOs! (1: 6),Wgts! (1: 6) 

1 

PRINT ’’FREE ARRAY & STACK AT STARTUP”, FRE(-l), FRE(-2) 
CALL Get Stat ions (A1 1st a^'K ) ,Cm0s! ( ) ,Wgts! ( ) ) 

Lprint CmOs! (1); CmOs! (2); CmOs! (3); CmOs! (4); ” CmOs” 

CALL SetFiles 

Print ’’Enter assumed position” 

Input ”AX=”,AP#(1) 
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Input "AY=",AP#(2) 

Input "AZ=",AP^'/(3) 

LPRINT "Assumed Position"; 

For i%=l to 3 

Print ^Assumed position in UTM (X,Y,Z)='* AP#(i%) 
LPRINT USING AP#(i%); 

Next i% 



‘j'c'jV’iV'vWc'sWc'sWrjV'jWc'sWoWcVoV'jWcVoWc'sWf'Wf^WvVc'WoWc'jV'WrVfVoV^WfyrVc* 



REM Aproximations for calling the subrut ine S0LVERANG2* 



SFAC!=1. 0 
SHGTy/=APy/(3) 

AP#(3)=0. do 

LHi'/=l. 1 
Rm#=6368459 



'Approximations for the subrutine for the 
' corrections that have already done 



'Convergence Limit test (meters) 
'Mean radius of Earth in meters 



* 0. 9996 



‘V:)VVrV?*‘V'VVc*'5V**>V'}VVfVcV?‘;VVr'5V'}VV?*VfVrVr‘V*VfVr‘jV*'5V* 



REM: Declaration of Weight matrix * 

' ^V*'5V^V*‘^VVrVvV-“V**Vc*Vc**V?V-^V*VrVfVf'5VVr*‘V'iVVr‘V*** 

For i%=l to N% 

W! (i%)=. 25 



'Diagonal matrix assuming 
' independent measurements 
Next i% 

' ‘,V^V****Vr^V’iV*Vr‘5VVr*'5VV<'V<’’:V**‘5V*:VV-‘V****'jV*‘/c’Vr*'5V 

Beep 2: CLS 

PRINT "Put Printer online in COMPRESSED PRINT Mode! " 

PRINT " Hit any key when ready to continue" 

WHILE NOT INSTAT 
LOCATE 2.37 
PRINT "_" 

WEND 

^ Print INKEY$ 

' LOOP TO READ AND PROCESS ALL DATA UNTIL END OF FILE HIT 

t 

'WHILE NOT EOF(l) 

For k°o=l to 500 

CALL READNICKMRC sec& , code%( ) ,R(/( ) ,X#( ) ,Y^f( ) ,ZK ) ,AIlSta#( ) ) 

REM Correction for slope distance * 

For i%=l to N% 

If (Ri'/Ci?^) > l.dO), Then 
DHGT#=Z(Ki%)-SHGT// 

HGT#=R(K i%) 2-DHGTy/ 2 

IF HGT,>‘< O.dO, Then 

PRINT "W-arning - at time sec&; " DHGT exceeds Range" 
PRINT " Hit any key to continue" 

While Not INSTAT 
Delay 0. 2 
Print 
Wend 
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t 



Print INKEY$ 

Else 

R//(i%)=SQR(HGT//) 
End If 

Else 



R#(i%) =0. dO 
End If 

Print "Ranges corrected for slope=" 
Next i% 



R#(i%) 



REM Calculation of the scale 



^ ix * •> V Vf Vc *iV 1 V Vr Vr * ' 



•jW? Vc* ^■jWr'jV * Vr ' 



factors * 

’ VoWr Vf * Vr * Vf* * 



AXP,-'/=AP,-'/(l) -500000 
For i%=l to N?ii 



LX#( i%)=X?'/( i?0 -500000 

M,-'Ki%)=AXPi'/ 2 + AXP^^lx/ZU^) + LX#(i%) 

SF//(i%)=0. 9996 + . 9996*M</(i%) /(6*Rm--2) 

Print "Scale Factors for UTM=" SF#(i/o) 

' Lprint ''Scale factors for UTM="SF/f ( i%) 

REM Correction for the scale factors * 



R//( i%)=R,^K i“0*SF//( i?o) 

Print "Corrected ranges=" R#(i?o) 
' Lprint "Corrected ranges=" 

Next i% 
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REM Formation of station position matrix 



For i/o=l to N% 
SP#(i?o,l)=X#(i%) 
P,^.‘(i?i,2)=Y;Ki%) 
SP;Ki%,3)=0. do 
Next i% 



CALL SolveRangX(AP//(),SP#(),R5'-‘(),Cm0s! (),N%,SFAC! , W! () , LM,^/ ,BRes#(),_ 
sxx! , syy! , sxy! , variance! ) 

CALL RiteSolCsec&jAPj-'rC ) ,code?o(3 ,Rj'/( ) ,BRes?A( ) ,sxx! >syy! ,sxy! .variance! ) 



LPRINT USING " 

For i% = 1 to 4 

LPRINT USING 

Next i% 

LPRINT USING " 

' WEND 
1 



nuim sec& AP?K1) AP-;'/(2); 

" y/##// #//#,-'/.# +M.//"; code%(i%) m±%) BRes#(i%); 
Mih inhyf. ir HHHK iHUHhU"', sxx! syy! sxy! variance! 



Next k?i 
CLOSE #1 
CLOSE //2 
END 

ERRORHANDLER: 

PRINT "FREE AT ERROR" ,FRE( -1) , FRE(-2) 

PRINT "AN ERROR TYPE", ERR, "HAS HAPPEN AT", ERADR 
END 
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* ***‘)V‘>V‘>VVf'>V‘jV'}V‘)V'>V‘)V->V*>V'>V')V‘)V'jV‘)VVr*-jVVr'»V'>V'>V'>V'>V'>V‘)V'>V‘)V'>V‘)V*')VVc’Vc'>V'>V‘5V'>V'>V^VVc‘>VVf*')V*Vc’Vr 

SINCLUDE "openfile. sub" 

'************* GetStations Subroutine ****Vf*Vc***** kjsiOct 89 * 
' file: NICKMR. SUB 

SUB GetStations(AllSta#(2),Cm0s! (l),Wgts! (1)) 

LOCAL i%, s$, code% 



CLS: PRINT "insert Disk with File ST. DAT in Drive B: 
Print " Hit any key to continue when ready"; 

While Not INSTAT 
Locate 2,37 
Print "! " 

Wend 



Print INKEY$ 

CALL 0PENFILE("I",1,"B: ST. DAT", 128) 

For i% = 1 to 6 

INPUT #1, code%,AllSta#(i”/i,l),AllSta#(i%,2),AllSta#(i%,3),Cm0s! (i%) ,_ 
Wgts! ( i?0 ,s$ 

Next i% 

CLOSE in 
END SUB 



************ SetFiles Subroutine 



*‘5V'5Wr7V')V*^VVr'>V9V'>WoWr'>WfVrVf kjsi Oct 89 ’’'Vc'jVVcVr'jVVc’ 
file: NICKMR. SUB 



^opens input data file of Miniranger values(as #1), and opens output file 
' for the computed results (as 



SUB SETFILES 

LOCAL s$ , txt$ ,FiName$ 

CLS: s$ = (Drive: path name, ext) 

PRINT "Enter MiniRanger data file for input" + s$; 

INPUT FiName$ 

CALL OPENFILE( "I" , 1 ,FiName$ , 128 ) 

PRINT "Enter solutions file for output" + s$; 

INPUT txt$ 

CALL 0PENFILE("0",2,txt$, 128) 

PRINT ir2 y "Solutions from file "; FiName$ 

END SUB 

^ ****vrvrvr**** ReadNickMR Subroutine ***vr**********%vv^^ kjs\ Oct 89 
’ file: NICKMR. SUB 

’specific to read mini-ranger files (from #1) from Krionertis thesis work 
’with time in secs & 4 station numbers & 4 ranges. 



SUB READNICKMR(sec&,code?o(l) ,rang//(l),X//(l),Y#(l),Z#(l),AllSta/i(2)) 
LOCAL s$ , i?i,nucode?o( ) 

DIM DYNAMIC nucode?o( 1: 4) 

I 

LINE INPUT in> s$ 
sec&=VAL(MID$(s$,8,7)) 

For i% = 1 to 4 

nucod e%( i°„ ) =VAL( M I D $ ( s $ , 1 1+4* i% , 4 ) ) 
rang#( i%)=VAL( MID$ ( s $ , 22+9*i% , 9 ) ) 
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'check nucode?o( ) , if not same as code%( ) then get new station coordinates 
I for the X#( ),¥,•'/() ,Z#( ) arrays 



If (code%(i%) <> nucode%( i%) ) , Then 
code%(i%) = nucode%(i%) 

X#(i%) = AllSta#(code%(i%),l) 

Yy/(i%) = AllStaf'/(code%(i%),2) 

Z#(i%) = AllSta#(code%(i%) ,3) 

End If 
Next i% 

END SUB 

' *-,v*********RiteSo I Sub rut ine **********v<-,v******' 



****■>•** kjs: Oct 89***** 
File: NICKMR. SUB 



Writes Miniranger Solution and error statistics to File y/2 



SUB RiteSol(sec6c,APy/( 1) ,code%( 1) 1) ,BRes?lf( 1) ,sxx! ,syy! ,sxy! ,vari! ) 

LOCAL i% 

I 

PRINT #2, USING "iH^rUU sec& AP#(1) AP#(2); 

f$ = "U mUihil +##.?'/" 

For i% = 1 to 4 

PRINT in, USING f$; code%(i%) R#(i%) BRes#(i%); 

Next 1% 

PRINT in, USING " M.m U-U +ifihU Hm.iiir-, sxx! syy! sxy! vari! 

END SUB 

'File = SOLRANGX. SUB 

>* * * * * * SOLVE RANGE POSITION (ver X) * * -> SUBROUTINE * * * 

I 

^Returns grid x,y in first 2 elements of apsnf/C ) 

^INPUTS: 

^apsn7'r(3) = approximate position as grid x,y,z 

^ spsns^r = station positions in 2D array as grid x,y,z 

' (nsta%,3) (array has nsta% rows by 3 columns for x,y,z) 

^ ORSfKnsta?^^ observed ranges to each station in spsns-:!^ row order 
’ rangcor! (nsta%) = C-0 corrections for each range reading 
*nsta% = number of stations/ranges 

^scalefac! = grid scale factor (grid; true ratio) 

* wgts! (nsta?o) = weights for each ran§^ (1/sigma 

’dxlim^'/ = convergence limit to stop iterating position 

^ solution (given in grid units - e. g. meters) 

^OUTPUTS: 

’apsn//(3) = x,y position returned in 1st two elements 

^ (elevation of antenna, z, unchanged) 

^ BRes 5 ^/(nsta%) = residuals on observed ranges ( ob-computed) 

*sxx! = variance in x (from Inv(ATWA)) 

^syy! = variance in y 

^sxy! = covariance of x & y 

Variance! = est of sigma-o (2?l^^g^:]B/4ii(i^-2) ) ) 

I 

$ INCLUDE "LUDCMP. BAS" 

$ INCLUDE "LUBKSB.BAS" 

SUB SolveRangX( apsn#( 1) ,spsns?!i(2) ,ORS,-'/( 1) ,rangcor! ( 1) ,nsta%,scalefac! , 
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wgts! ( 1) jdxlimj'/jBResj'K 1) >sxx! ,syy! ,sxy! , variance! ) 
LOCAL A^K ) , spsn/K ) ,ORa# , Ax# , Ay# , Bi# , ATW#( ) , ATWB#( ) , ATWA#( ) , Indx%( ) , 
sum#, iteration%, i%, j%,k%,YCol#( ) , ATWAInv#( ) ,xdx%( ) , Btemp#( ) 
DIM spsn#(l:3), ATWB#(1:2), ATWA#( 1; 2 , 1: 2) , Indx%( 1: 2) ,xdx%( 1: 4) 
Btemp#( 1: 4) 

T — 

i%=l 

For j% = 1 to nsta% 

If ORs//(j%) > l.dO, then 

xdx%(i%) = j% *xref to non-zero ranges (range > Im) 

INCR i% 

End If 
Next j% 

nranges% = i%-l 

DIM A//( 1: nranges%, 1: 2) , ATW//( 1: 2 , 1: nranges%) 

SolveAgain: 

For i% = 1 to nranges% 

ORa//=ORS//(xdx%(i%)) + rangcor! (xdx%( i%) ) 

For j/o = 1 to 3 

spsn7Kj%) = spsns//(xdx%( i%) , j%) 

Next j% 

CALL RangeLOPC apsn//( ) ,spsn#( ) jscalefac! ,ORa//,Ax//,Ay//,Bi//) 

AtK i%, l)=Ax-:'/ 

A7'/(i%,2)=Ay7^.^ 

BRes/r( i%)=Bi// ^Residual (o-c) on each range 

Next i% 

I 

’ Compute ATW 02 ,nranges?o | , & ATWB 02 | 

For j% = 1 to 2 
sum//=0. dO 

For i% = 1 to nranges/b 

ATC//(j%,i°o) = A//(i?o,j%)*wgts! (xdx%(i?0) 
sum?'/ = ATWf/( i%) -'BRes7*'/( i?o)+sum7i^ 

Next i% 

AT\vB//( j%)=sum 7 '/ 

Next j?o 

t 

* Compute ATWA 02,2 | 

1 

For j/o = 1 to 2 
For i% = 1 to 2 
sum7’^r = 0. do 

For k% = 1 to nranges% 

sum^/ = ATVvV/( j%,k%)*A//(k%,i%)+sum// 

Next k% 

ATCAy/C j%,i%)=sum// 

^Next 

CALL LUDCMPC ATCA//C ) , 2 , Indx%( ) , i%) 

CALL LUBKSBC ATWA//C ) , 2 , Indx%( ) , ATOB//C ) ) 

sum 7 *!-= 0 . do 

For i®b =1 to 2 

apsnrrC i%)=apsnf(K i%)+AT\vB//( i%) 'shifts position by dx,dy which has 
sum// = ATWB?2( i*°o^um# 'been returned as ATWB() from LUEKSB 
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Next i% 

sum^/ = SQR(sum//) 

If sum// > dxlim//, Then 
INCR iterationi 
If iteration% > 7, Then 
Beep 3 

Print ^^Position did not converge --NO SOLUTION! 

Print Using ’^Delta Position = +//. ////////meters**; sum// 
sxx! = 0.0: syy] = 0.0: sxy! = 0.0: var! = L0G10(sum//) 

EXIT SUB 
End If 

GoTo SolveAgain 
End If 

t 

* Compute variance-covariance values from inverse of ATWA 

* sxx! = sigma-x 2 

* syy! = sigma-y 2 

* sxy! = sigma-xy 

ERASE ATW//, spsn// *make some room in memory 

DIM YCol//(l:2), AT\v'AInv//( 1: 2 , 1: 2) 

For j% = 1 to 2 

For i?^ = 1 to 2 * identity vector, YCol, to find 

YCol 7 '/(i%) = 0. do 'inverse using back-substitution 

Next i% 

YCol//( j°o) = l. dO 

CALL LUBKSBC ATWA//C ) , 2 , Indx?o( ) , YCol//( ) ) 

For i?o =1 to 2 

ATl\’A I n v// ( i% , j % ) =YC o 1// ( i% ) 

Next i% 

Next j% 

sxx! = CSNG(AT\\^AInv//(l,l)) 
syy! = CSNG(AT\vAInv/^(2,2)) 
sxy! = CSNG(ATWAInv//(l,2)) 

I 

* Compute estimate of unit 

* Note that range residuals are in BResiK) 

I 

If nranges% = 2, Then 

variance! = 9999 'signifies no variance computable 

Else 

sum// =0. do 

For i% = 1 to nranges?o 

sum// = Bres//(i%) 2-wgts! (xdx%( i%) ) + sum/^ 

Next i% 

variance! = CSNG( sum///(nranges%-2) ) 

End If 

» 

If nranges?o < nsta°o, Then 'Set residuals into proper channel 

For i?o = 1 to nsta®o 
Btemp//(i%) = 0. dO 
Next i% 

For j% = 1 to nranges% 

Btemp?/(xdx/o( j%) ) = Bres//(j%) 

Next j°o 

For i% = 1 to nsta% 
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If nranges% = 2, Then Bres#(i%) = 0. dO 
Bres#(i%) = BtemptfCi^) 

Next i% 

End If 

I 

END SUB 

I 

Sv * * * Vf * RANGELOP * * SUBROUTINE * * * * * * * * * 

I 

* Computes A matrix components Ax & Ay and grid range 
^ for input to multiple-LOP solution algorithm. 

* Inputs: apsn#() = approx position x,y,z coordinates (grid) 

' spsnf'/C) = station x,y,z coordinates 

* scalefac! = grid scale factor (point or line) 

' ORa// = observed slant range 

^Output: Ax# 6c Ay# such that (Ax)dx+( Ay)dy = OR-CR 
' OmCR,# = Observed minus Computed range at apsn# 

^Note: Observed range (OR) has been corrected for elevation 
' differences and grid scale before comparing with CR. 



SUB RangeL0P(apsn#( 1) ,spsn#( 1) jScalefac! 

ORa# , Ax# , Av# , OmCR# ) 

LOCAL i%,CR#,dx#() 

DIM dx#(l: 3) 

I 



For i% = 1 to 2 

dx#(i%) = apsn#( i?i) -spsn,^/( i/4) 

Next i% 

CR,v‘ = SQR(dx,-'Kl) 2+dx#(2) 2) 

Ax# = dx#( l)/CR,-;i 
Ay# = dx#(2)/CR# 
dx#(3)=0. do 

OmCR# = ORa#*scalefac! - CR# 

END SUB 

k *V V: Vr Vr Vr Vr Vc V? • V Vr :V * Vr Vc Vr Vr ■jV Vf Vc V? Vr * Vc Vf Vr V? * V: Vf V? V? Vc Vc Vc Vr Vc * Vc Vc Vr Vc V? * -V V? V? * Vr Vc V? Vr V? -V -V V: Vr V- V? V? Vr Vr Vr V: Vr 
V Vr Vr V? Vr Vc V? V? Vr «V Vr ‘V Vr Vr V? ‘V Vr ‘V -V Vr *V nk Vr *V *V *V Vf Vr Vr Vc Vr *V V r Vr V? Vc Vr *V * V Vr */c Vr y.' *V Vr Vr Vr *A* nk Vr Vc - c *V Vr *A* Vr ‘V 



file: OPENFILE. SUB 



Opens a disk file using TBasic commands. Preforms error checking and 
allows recovery from common mistakes. 

INPUTS: modes$ = One character designation for filemode 

”0^' = for sequential output (to be written to) 

”l” = for sequential input (to be read from) 

Csee Turbo Basic Ref, pg 282, for other choices | 

FileNum% = file number to be used in read/write calls 
FileName$ = Drive: /path/name. ext of file 
RecLength?b = record size in bytes, 1 to 32767, typical = 128 



SUB OpenFile(mode$ ,FileNum?o,FileName$ ,RecLength%) 



LOCAL errnum°o, ans$ 
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OpenSequence: 

ON ERROR GoTo ErrorHandling 

OPEN mode$ ,FileNum/o ,Fi leName^ ,RecLength% 

ON ERROR GoTo 0 
EXIT SUB 

f 

ErrorHandling: 

Beep 2 

errnum% = ERR 
Select Case errnum% 

Case 53,75,76 

Print "FILE or DRIVE PATH Not Found for . , . " FiIeName$ 

Input "Re-Enter Drive: path name, ext > "; FiIeName$ 

Case 68,71 

Print "DISK DRIVE NOT READY for ", FileName$ 

Case Else 

Print "Opening File Gave Error Number ", errnum% 

End Select 

Input "(R)etry or (H)alt Program "; ans$ 

If Ucase$(ans$) = "H" , Then 
Resume HaltSequence 
Else 

Resume OpenSequence 
End If 



HaltSequence: 



'terminates main program 



Print "Program Halted! ! ! " 

Print "Subroutine OpenFile with Error Number ",errnum% 
Print " encountered at program address ", ERADR 
END 

END SUB 









' LU Back Substitution 

I 

' Ref: Numerical Recipies, by WH Press, page 37 
\ 

SUB LUBKSB ( A,"/ ( 2 ) , N?o , INDX%( 1 ) , B#( 1 ) ) 

LOCAL I I?i yi7o, j%, sum# 

I 

ii%=0 

For i%=l to N% 

II%=Indx%( i%) 

sum#=B#(II?0 

B#(II%)=B#(i%) 

If ii%<>0, then 

For j%=ii% to i%-l 

sum#=sum#-A#( i%, j%) 

Next j?o 

Elself sum#<>0, dO Then 



End If 
B#( i°b)=sum# 



Next i% 



' ****v 

' I'rMri'ti 



For i%=N% to 1 Step -1 
sumi^=BiKi"/o) 

If i%<N%, Then 

For j%=i%+l to N% 

sum#=sum5'/-A#( i%, j%) 

Next j/o 
End If 

Next i% 

END SUB 

'c^VVc*VcVf*yr*VrVr**V-‘>VVfVr*‘jV';V‘sV**-5V'>V->V****->V'>VVr*>V->V'5V'>VV?‘>VVf'jV****VfVf‘>V*Vr*^VV< 

V*^VV?Vr*‘5VVf*Vr*Vr**Vr^V*‘jVVr->VVr**'jV'5VVfVc')Vyc'5V':V***5VVf-5V**'sV*'5V*/V‘;V***VryfV?*VcV< 

’ LU Decomposition Subroutine 

^ Ref: Numerical Receipes, by WH Press, pages 34-36 



*Vr*" 



****** ^’ 
V***VoW 



SUB LUDCMP(A//(2) ,N%,Indx?a 1) ,D%) 

LOCAL vv//( ) ,Dum 7 *^/, i%, j%,k/o, imax%,aamax//,sum# 
STATIC nmax?i, tiny?"/ 



nmax%=10 
tiny//=l. OD-20 



DIM vVt'/C 1: nmax?i) 



D%=1 

For i%=l to 
aamax//=0 
For j%=l to N% 

I f ABS( At'/C i°o, j%) )>aamax// , Then aamax//=ABS( A#( i?b, j%) ) 
Next j/o 

If aamax//=0, Then 

Print ’’Matrix has zero row // ” i% 

Exit Sub 
End If 

vv?Ki%)= 1/aamax# ’save row scaling factor 

Next i/o 

t 



For j%=I to N% 

If j%>l, Then 
For i/o=I to j%-l 

Sum^/=A^/( i/o, j°i) 

If iVl Then 
For k?i=l to i%-l 
Sum-^‘=Sum;'/-A,-',f( i%,k%)*A/Kk%, j%) 

Next k% 

A#( i%, j%)=Sum,''/ 

End If 
Next i% 

End If 

' continuing j loop, searching for largest pivot 

aamaxri =0 

For i%= j% to N°o 
Sum;/=A,"'/( i%, j?i) 
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If j%>l, Then 
For k%=l to j%-l 
Sum#=Sum#-A#( i%, 
Next k% 

End If 

Dum#=vv,''/( i%)*ABS( Sum,-^/) 
If Dumj'/>aamax,'f , Then 
imax%=i% 
aamax#=Dunn'/ 

End If 
Next i% 



*A#(k%,j%) 



'figure of merit for pivot 



If j%<>imax%, Then 
For k%=l to N% 

Dum#=A#( imax%jk%) 

A// ( imax% , k% ) ( j % , k% ) 

A#(j%,k%) = Dum// 

Next k% 

D%=-D% 

vviK imax%)=vv#( j%) 

End If 



Indx%( j%)=imax°o 
If (j%oN%), Then 

If AH}%,r/o)= 0, Then A#( j%, j%)=tiny# 
Dum,-'/= 1/MK 3%, 3%) 

For i%=j%+l to N% 

m i%, j%)=A,-'K i%, j%)*Dum,-'/ 

Next i% 

End If 
J » 



'interchange rows if reqd 



' change parity of D% 



'divide by pivot 



If A#(N%,N%)=0, Then A#(N?o,N?4)=tiny# 
END SUB 



* **■! 
1 1 t 


^cVr**-V***VrVrVr.W 


'.-******-.>****** 


VfVrVrVf’jWrV 


r***Vr, 






sWr Vr VrVfVr Vr *V' 


1, 


585260. 161, 


4092490. 284, 


51. 7, 


0. 0, 


0. 


25, 


TREVOR 


2, 


602947. 672, 


4085231. 399, 


23. 1, 


0. 0, 


0. 


25, 


WATS 


3, 


609863. 128, 


4076611. 345, 


33. 9, 


0. 0, 


0. 


25, 


PACKARD 


4, 


607621. 289, 


4055915. 264, 


137. 2, 


0. 0, 


0. 


25, 


HAYES 


5, 


600434. 009, 


4051260. 014, 


10. 0, 


0. 0, 


0. 


25, 


DOPPLER 


6, 


596669. 490, 


4051826. 411, 


134. 5, 


0. 0, 


0. 


25, 


HANK 
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APPENDIX B. GEOMETRICAL ACCURACIES OF THE MONI EREV 

BAV MINI-RANGER NETWORK 




Figure 13. SIAllON COMBINATION TREVOR- WAT S-IIAVS-I lANK 
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STATION COMBINATION 

TREVOR-PACKARD-DOPBLER-IIANK 
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figure 15. STAIION COMBINATION WAIS-PACKARD-IIAVS-HANK 
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Figure 16 . STAIION COMBINATION PACKARD-MAVS-DOITLEIMIANK 
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Figure 17. STATION COMBINATION TRE\ OR- WATS-IIANK 
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Figure 18. STATION COMBINATION T REVOR-PACKARD-I lANK 
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Figure 19. STATION COMBINATION WATS-FIAVS-HANK 
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Figure 20. STATION COMBINATION TREVOR-PACKARD-HAVS 
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